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Abstract 

We solve the inverse scattering problem for multidimensional vec- 
tor fields and we use this result to construct the formal solution of 
K* . the Cauchy problem for the second heavenly equation of Plebanski, 

^1 I a scalar nonlinear partial differential equation in four dimensions rel- 

^^ ■ evant in General Relativity, which arises from the commutation of 

'^ . multidimensional Hamiltonian vector fields. 

O 

\o 

O 

1 Introduction 

Ch ■ In this paper we solve the inverse scattering problem for multidimensional 

K^ ■ vector fields and we use this result to construct the formal solution of the 

K^ , Cauchy problem for the real second heavenly equation 

^ 0t.-0,y + e,Jyy-9ly = 0, 9 = 9 {x , v , z , t) E R , x,y,z,teR, (1) 

where subscripts denote partial derivatives. 

This scalar nonlinear partial differential equation (PDE) in 4 indepen- 
dent variables x, y, z, t, introduced in pi by Plebanski, describes the Einstein 
field equations that govern self-dual gravitational fields. Its 2-dimensional 
reduction 6'^ = ^t = is the Monge- Ampere equation, relevant in Differential 
Geometry. As we shall see in the following, the heavenly equation plays also 



a distinguished role in the theory of commuting, 2-dimensional, Hamiltonian 
dynamical systems. 

The heavenly equation, together with the equations for the self-dual Yang- 
Mills (SDYM) fields 0, are perhaps the most distinguished examples of non- 
linear PDEs in more than three independent variables arising as commuta- 
tivity conditions of linear operators [3], and therefore amenable, in principle, 
to exact treatments based on the spectral theory of those operators |ll,j3]. 
If the SDYM equations are considered on an abstract Lie algebra, then the 
heavenly equation can actually be interpreted as a distinguished realization 
of the SDYM equations, corresponding to the Lie algebra of divergence free 
vector fields independent of the SDYM coordinates jn|. Equation (jT)) has 
been investigated within the twistor approach in [7j,[H]. A bi-Hamiltonian 
formulation and a hodograph transformation for (^ have been recently con- 
structed in [£j and ^10, ; a nonhnear 9-dressing and a generating equation for 
its hierarchy can be found in [TT] . 

The present paper is a simplified version of the manuscript 



2 Integrable PDEs in arbitrary dimensions 
and the heavenly equation 

It is known that the commutation of multidimensional vector fields leads to 
nonlinear first order multidimensional PDEs (see, f.i., [3]). In this spirit, we 
derive now a basic class of integrable nonlinear PDEs in arbitrary dimensions 
possessing, as distinguished reduction, the heavenly equation (jT)). 
Consider the following pair of vector fields 

N 

L, = dt, + \d,^ + Y, ^%, = 9u + A5., + u,-V^, z = 1, 2 (2) 
fc=i 

where dx denotes partial differentiation with respect to the generic variable 
X, X = {xi,..,xn), V^ = {dx^,-.,d^j^), Ui = {u},..,uf^), i = 1,2, A is a 
complex parameter and the vector coefficients Ui depend on the independent 
variables ti,Zi,Xk, i = 1,2, k = 1,..,N, but not on A. The existence of a 
common eigenfunction / for the operators Li and L2: 

Lif = Uf = 0, (3) 



implies their commutation, VA: 

[^1,^2] = 0, (4) 

which is equivalent to the following determined system of 2N first order 
quasi-linear PDEs in (4 + A^) dimensions for the 2N fields •ui,U2'- 

Ult2 - ^"2ii + {U2 ■ Vs) Ul - {Ul ■ Vs) M2 = 0. 

Parametrizing the first set of equations in terms of the vector potential U 

Ui = Uz,, i = l,2, (6) 

one obtains the following determined system of A^ nonlinear PDEs for the A^ 
dependent variables [/ in (4 + A^) dimensions: 

Ut,z2 - Ut2z, + {U., ■ V,^) U,, - {U,, ■ V,) tJ,, = 0. (7) 

This system admits a natural reduction; indeed, applying the operator 
V^- to equations ((Tj), one obtains 

dt,d,, - dt.A, + (Uz, ■ v.^) d,, - [tj,, ■ v.^) a,,] (v.. ■ t/) = 0, (8) 

from which one infers that the condition 

Vs-tJ = (9) 

is an admissible reduction for equation ((Tj), implying that the condition of 
zero-divergence : 

Vs-Ui = 0, 2 = 1,2 (10) 

is an admissible constraint for the vectors -Uj, i = 1,2. 

From now on, we concentrate our attention on the following important 
example: 

AT = 2, Zi = Xi, i = l,2 (11) 

and we make the following change of notation for the remaining 4 variables: 

ti = z, t2 = t, xi = X, X2 = y. (12) 



Then the system (jTj) reduces to the following determined system of two PDEs 
in 4 dimensions: 



Utx - Uzy + (Uy ■ V,.) U, - (U, ■ V,.) Uy = 0, 

UeR^ x={x,y), Vs=id,,dy), 



(13) 



corresponding to the Lax pair: 



Li = d^ + \dx + ui ■ Vg, Ml = U^, 

L2 = dt + Xdy + U2 ■ Va;, U2 = Uy 



(14) 



In this case, the zero-divergence reduction ^T^ makes the two vector fields 
Ui-Vsc Hamiltonian, allowing for the introduction of two Hamiltonians Hi, i = 
1, 2 such that: 

Ui = (Hiy, -Hi^), z = 1, 2 (15) 

which, due to (jS^), are parametrized by a single potential 6: 

H-i = Uxi, U = [Oy, —Ox), /-jg\ 

Ul = [Oxyi ~Oxx)i ^2 = [vyy, —Oxyj- 

Then the compatible linear problems ©,(^3) can be written down as Hamil- 
ton equations with respect to the times z, t: 



ft = {H2-Xx,f},, ^"' 



where {■, -j^ is the Poisson bracket with respect to the variables x, y: 

{f,9}£ = fxQy- fyQx, (18) 

and the nonlinear system (jl3|) reduces to the heavenly equation in Hamilto- 
nian form 

Otx - Ozy + {Ox, Oy}£ = constant, (19) 

equivalent to (jT)) after choosing the constant to be zero. 

3 1ST for the nonlinear PDEs dH) and (^ 

Since the Lax pair (fT^ is made of vector fields, Hamiltonian in the heavenly 
reduction (fTUI) . (fTTj) . the eigenfunctions satisfy the following basic properties. 



which will introduce important novelties in the Inverse Scattering Transform 

(1ST). 

1) The space of eigenf unctions is a ring: if /i, /2 are two solutions of the 
Lax pair ()14|) . then an arbitrary differentiable function F(/i, /2) of them is 
a solution of (J14p . 

2) In the heavenly (Hamiltonian) reduction mU\} . ill]) , the space of eigenfunc- 
tions is also a Lie algebra, whose Lie bracket is the natural Poisson bracket 
M8\j : if /i, /2 are two solutions of the Lax pair (jl4|) . then their Poisson 
bracket {/i, f2\{x,y) is also a solution of (HD). 

Now we consider the Cauchy problem for the system ()13|) and for its 
heavenly reduction ((T]) , within the class of rapidly decreasing real potentials 



ui: 



u\ -^ 0, (x2 + 2/2 + z^) -^ oo, . . 

u\^^, (x,y,^)GM^ t>0, ^ ^ 

interpreting t as time and the other three variables x, ?/, z as space variables. 
To solve such a Cauchy problem by the 1ST method (see, f.i., |4j,|5j), we 
construct the 1ST for the operator Lx in ((H^), within the class of rapidly 
decreasing real potentials, interpreting the operator L^ in ((TUd) as the time 
operator. The formalism presented here can be generalized to the vector 
fields (121) in a straightforward way. 

3.1 Basic eigenf unctions 

The localization (f^UI) of the vector potential u\ implies that, if / is a solution 
of L\f = 0, then 

/(f,2;, A) -* /±(^, A), 2;^±oo, .^x) 

f :=f- (A, 0)2;; 

i.e., asymptotically, / is an arbitrary function of (x — Xz), y and A. 

A central role in the theory is played by the real Jost eigenf unctions 
(pi,2{x, z, A), the solutions of I/iv?i,2 = uniquely defined by the asymptotics 

ipi{x,z,X) ^ X - Xz = ^, ip2{x,z,X) ^ y, z ^ -00. (22) 

In this paper we often use the compact vector notation: / = (/i, f2)^- Then: 



The Jost eigenvector is equivalently characterized by the integral equation 

/jg3 dx'dz'G{x — x', z — z'; A) ('Ui(x', z') ■ Wg') 0{x', z', A) = ^, 

in terms of the Jost Green's function 

G{x,z;X) = e{z)6{x-Xz)5{y). (25) 

A crucial role in the 1ST for the vector field Li is also played by the 
analytic eigenfunctions V^±(x, z, A), the solutions of Li'i/j± = satisfying the 
integral equations 

ij±{x, z, A) + /]g3 dx'dz'G±{x — x' , z — z'\ A) {u\{x' , z') ■ V^) ip±{x', z', A) = ^, 

(26) 
where G± are the analytic Green's functions 

The analyticity properties of G±{x, z, A) in the complex A - plane imply that 
il^+{x, z, A) and i/j-{x, z, A) are analytic, respectively, in the upper and lower 
halves of the complex A - plane, with the following asymptotics, for large A: 

4(f,^,A)=e+^±^ + 0(A-2). |A|»1, (28) 

A 

where: 

(X oo\ 

l-l\dx'u,{x',y,z), 
-OO X / 

(29) 
entailing that 

Ml (f , z) = -Q±x (f, z) . (30) 

It is important to remark that the analytic Green's functions (J27|l exhibit 
the following asymptotics for z -^ ±oo: 



G'±(f,2;A)^±2^^g-j, z^+oo, 
G±(f,2;A)^±2;^g-i, z^-oo, 






It follows that the z = +cxd asymptotics of i/j^ and ip- are analytic respec- 
tively in the lower and upper halves of the complex plane ^, while the z = — cxo 
asymptotics of -0+ and ip_ are analytic respectively in the upper and lower 
halves of the complex plane ^. This mechanism was first observed in jTB]. 



3.2 Spectral data 

The z = +CXD limit of tp defines the scattering vector a of Li: 

hm 0{x, z- A) = S{1 A) = i + a{l A) . (32) 

The direct problem is the mapping from the real vector potential mi, 
function of the three real variables (x, 2), to the real scattering vector a, 
function of the three real variables (^, A). Then the counting is consistent. 
The impact of the Heavenly constraint (jl(Jj) . (jlip on the spectral data will be 
discussed in Section 3.5. 

The Jost solutions </?i,2 form, together with the constant eigenfunction 
A, a basis in the space of the eigenfunctions of Li (which is a ring). The 
representation of the analytic eigenfunctions 'ip± in terms of (p yields: 

^±(f, z, A) = /C± ((^(f, z, A), A) = 0{x, z, A) + x± (<^(^, z, A), A) , (33) 

and this formula defines the spectral data x±- Since the z -^ —00 limit of 

dSni) reads: 

hm 4-? = X±(?'',A), (34) 

2:— >— 00 

the above analyticity properties of the LHS of ()34j) in the complex ^ - plane 
imply that x+(^) A) ^^^ X-(?! ^) are analytic respectively in the upper and 
lower halves of the complex plane ^, decaying at ^ ~ 00 like 0{^~^). There- 
fore their Fourier transforms x+i'^i -^) a-^d X-{^i ^) have support respectively 
on the positive and negative real ui semi-axes. 

The spectral data x± can be constructed from the scattering vector a 
through the solution of the following linear integral equations 



(35) 



X+{uJ, A) + e{ui) \^a{uj, A) + /jg2 df] x+iv, >^)Qiv, uj, A) j = 0, 
X-{uJ, A) + 9{-u}i) [a{u}, A) + /jj2 df] x-{v, X)Q{f], u), A) j = 0, 

involving the Fourier transforms a and x± of c? and x±'- 

~a{CJ,X)= [ d^ai^,X)e-'^<, ^±(^,A)= / d^x±€ X)e-'^< (36) 

Jr2 Jr2 



and the kernel: 



Q{fj, u,X)= I ^e'(^^--^)-?V'^-^(r.A) _ 1]. (37) 
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To prove this result, one first evaluates (|33|) at z = +00, obtaining 

lim 4 - e = <?(e, A) + x±(e + ^{^, A), A). (38) 

2:^00 

Applying the integral operator J^2 (i^e~*'^'>- for u;i > and oji < respec- 
tively to equations ((211)+ and (jHHjl -. using the above analyticity properties 
and the Fourier representations of x± and a, one obtains equations (jH^j) . 
We end this section remarking that the reality of the potentials: ui G M^ 

implies that, for A G M, (^= 0, V'+ = "0-! consequently: a = a, x+ = X-- 

3.3 Inverse Problem 

An inverse problem can be constructed from equations ()33|) . Subtracting ^ 
from equations (jHSJ- and ((221)+, applying respectively the analyticity pro- 
jectors J+ and P_: 

^ _ 1 /■ dA' 

^^ = ^2^y^A'-(A±.6)- ^^^^ 

and adding up the resulting equations, one obtains the following nonlinear 
integral equation for the Jost eigenfunction 0: 

(^(f, z,X) + ^J^ v^fri^X- i0i^^ 2;, A') , A') - 
2^/Mv3t^x4^1x,^,A'),A')=e ^ ^ 

Given the spectral data x±) one reconstructs the eigenfunction tp from ()40|). 
the analytic eigenfunctions from ()33p. and {ti from equation ()30|) . This in- 
version procedure was first introduced in ^^. 

3.4 ^-evolution of the spectral data 

As the potentials ui^2 evolve in time according to equation (fT^ . the t- 
dependence of the spectral data a and x±) defined in (fH^ and (jHSl), is de- 
scribed by the equation: 

a(e,A,t) = a(e-(0,A)t,A,0), .^^. 

X±(f,A,t) = x±(?"'-(0,A)t,A,0). 

To prove it, we first observe that 

(t)i{x,z,X,t) = ipi{x,z,X,t), (t>2{x,z,X,t) = (p2{x,z,X,t) - Xt (42) 

8 



are, together with A, a basis of common Jost eigenfunctions of Li and L2. 
The y = +00 hmit of equation I/2(;/> = yields at + XBy = 0, whose solution 
is dH^). Analogously, 

7T±^{x,Z, X,t) = iIj±i{x, Z, X,t), 7r±2{x, Z, X,t) = ij)±2{^i ^1 ^jt) ~ '^t (43) 

are a basis of common analytic eigenfunctions of Li and L2; therefore 

7r±i =/C±i((/), A), vr±2 = /C±2(0, A), (44) 

for some functions /C±i^2 depending on (x, z, t) only through (p. Comparing, 
at t = 0, these equations with equations (jH^ . one expresses /C-i-1^2 in terms 
of /C±^ 2 5 obtaining equations pib). 

3.5 The heavenly reduction 

In the heavenly (Hamiltonian) reduction ()10|) . ()11|) . the transformations ^ — >■ 
^{C: ^)) ? ^ •'^±(^5 ^) are constrained to be canonical: 

{5i,52}^-={/C±i,/C±2}g = l, (45) 

or, in terms of a(^, y. A) and x±(C, ?/, A): 

f^l? + f^2y + {f^l, ^^2}^- = X±^ + X±2y + {X±l, ^±2}^ = 0" (46) 

To prove it, we recall that the Poisson bracket of the eigenfunctions Lpi and 
ip2 is also an eigenfunction: Li{ipi,Lp2\x = 0. Using the asymptotics (j221), 
one infers that {'^i,'^2}x —^1, at z -^ — oo; therefore, by uniqueness, 
{v^i)¥'2}x = 1- Evaluating now this Poisson bracket at z = +cxo and using 
(|H^. one obtains the constraint (P3|) for S. We also observe that the eigen- 
functions {'?/'+ 1, '0+2 }x and {'?/'_ ]^, ■0_ 2 }x are analytic in the upper and lower 
halves of the A-plane and go to 1 at |A| ^ cxd. Since 1 is also an eigenfunc- 
tion, by uniqueness they are identically 1: {'4'±ii'4^±2}x = 1- Therefore, from 
the equations: 

{^±1, '0±2}x = {/C±i, IC±2}{v>i,v'2)Wu 'f2}x = 1, (47) 

consequence of (jH^ . one infers the constraints (jI3|) for /C±. 
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4 Commuting A- families of dynamical systems 

It is well-known (see, f.i., PSj) that linear first order PDEs like (jHI) are inti- 
mately related to systems of ordinary differential equations describing their 
characteristic curves. The vector fields Li,2 (fT^ are associated with the 
following two A-families of commuting dynamical systems: 



Li: f = (A,0)+Mi(f,z,t), 



L2: f = (0,A)+M2(f,^,t). 
In the heavenly reduction (fTUI) . (fTT|) . they read: 

Li: f + {i/i + A?/,f},^ = 0, 



(48) 



(49) 



L2: f + {//2-Aa;,f},^ = 0. 



Therefore: the heavenly equation characterizes the commutation, VA, of two 
X-families of Hamiltonian dynamical systems, with Hamiltonians: 

Hi = Hi{x, z, t) + Xy, H2 = H2{x, z, t) - Xx. (50) 

There is a deep connection between the above 1ST and the z-scattering 
theory for the commuting fiows (PHj) and ()49|) . Let (p{x, z, X, t) be the common 
eigenfunctions of Li and L2 defined in (H^ : then, solving the system uj = 
(f){x, z, A, t) with respect to x and y (assuming local invertibility) , one obtains 
the following common solution of the commuting fiows (j48|) : 

u; = (/)(x, 2;, A, t) -v^ x = f{z,t,X,u) ~ I , I + cJ, 2; ~ — 00. (51) 

The z = +00 limit of the solution f{z, t, A, u): 

x = r{z,t,X,u) ~ I . \+(1{lu,X), y ~ -|-oo (52) 

defines the scattering vector A{u) = Q{uj) —ujoi PHj) . which is connected to 
the 1ST data S by inverting the system u = S{x — A^, y — At, A) with respect 
to X and y: 

u = S{x-Xz,y-Xt,X,0) ^ x= ( ^^ j +Q{u,X). (53) 

10 



In the heavenly reduction, the transformation uj -^ f2(d;, A) is clearly canon- 
ical: {r2i,i72}(j = 1, and the scattering vector A exhibits the following con- 
straint: 

Ai^^+A2^, + {Ai,A2}^ = 0. (54) 

The above 1ST theory allows one to reconstruct, from the scattering vector 
A, the potentials u\^2 of the dynamical systems (jlHj) and (PUJ) . 



5 Other inverse problems 

Due to the ring property of the space of eigenfunctions, there are other 
distinguished ways to do the inverse problem. 

5.1 A nonlinear RH problem 

We begin with a more traditional (nonlinear) Riemann-Hilbert (RH) prob- 
lem. Solving the algebraic system (jHH|) _ with respect to Cp: = L{tlj_,X) 
(assuming local invertibility) and replacing this expression in the algebraic 
system ((221)+, one obtains the representation of the analytic eigenfunction 

— * — * 

"0+ in terms of the analytic eigenfunction -i/).: 

^+ = :^(^_,A) =^_ + ^(^_,A), AgM, (55) 

which defines a vector nonlinear RH problem on the real A axis. The RH data 
R are therefore constructed from the data x± by algebraic manipulation. 
Viceversa, given the RH data R, one constructs the solutions 'ip± of the 
nonlinear RH problem (j^3jl and, via the asymptotics (j^H|) . the potential ui. 
As for the other spectral data, one can show that the t-dependence of R is 
described by i?(f , A, t) = R{^ — (0, A)t, A, 0), and the reality constraint takes 

the following form, for A G M: 7?(7?(C, A), A) = C, Vf. At last, the heavenly 
constraint reads {T^i, 7^2}^ = 1, or, in terms of R{(, A): 

Ria + R2c, + {Ri,R2}^=0. (56) 

5.2 Linearization of the inverse problem 

It is possible to construct a linear version of the inverse problem of Section 3.3 
by exponentiating the Jost and analytic eigenfunctions used so far. Consider 

11 



the following scalar functions: 

<l>(f , z, A, a) = e^°-^(^^'"'^), ^±(f, z, A, a) = e*"-^^±(^^'"'^\ a G M^. (57) 

Due to the ring property of the space of eigenfunctions, also $ and ^± are 
eigenfunctions; $ is characterized by the asymptotics $ -^ exp{id ■ ^), z ^ 
— cx), while \E'±(x, z, Aa) are analytic respectively in the upper and lower 
halves of the A plane, with asymptotics: 

v^i = e-^-«^ (l + 'll^j^M + 0{X-^)] . (58) 

From equations ^7\i and (jH^ it follows that 
^±(f , z, A, a) = $(x, 2, A, a) + e^'^-'^^^'^'^) (e^-x±('?(--,^,A),A) _ ^^ _ ^^g^ 

The last term of this equality is clearly equal to its anti-Fourier transform: 

^ia-0 (^^id-x±{^,x) - 1) = /" d(3K±iaJ,X)e'^-^, ^ e ^\ (60) 

where 

K^{a, /5, A) ^ 4, ^e^(-/3>e[e-x±«.A) _ i]. (gi) 

Then, using (|37k). the last term in the RHS of (j3^ becomes a linear ex- 
pression in terms of $, and equations (jKn|l give the Hn ear representations the 
analytic eigenfunctions \E'± in terms of the Jost eigenfunction $: 

'$±{x,z,X,d) = ^{x,z,X,a) + d(3K±{d,(3,X)^{x,z,X,P). (62) 

Jr2 

Multiplying the equations (jU^ + and (j^ _ by exp{—id ■ ^), subtracting 1, 
applying respectively P_ and P+, and adding the resulting equations, one 
obtains the following linear integral equation for $: 

^(^^ «) + 2^ /m ITzfn^ L^ dpK.ia, /3, A')$(A', /3)e-^(^'-^)^- 

(63) 
- 2^ /m v^t^ /m^ ^^X(5, P, X'MX', /3)e-(^'-^)^ = e-^-^ 

in which we have omitted, for simplicity, the parametric dependence on [x, z) 
of $. Once $ is reconstructed from (j63|) and, via (j62p . \E'-|- are also known, 
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the potentials are reconstructed in the usual way from the asymptotics (|^ 
of ^±. 

We finally observe that the reality constraints for the eigenfunctions 
$, \E'± and for the data K± read, for A G M, 



^(x, z, A, a) = $(x, 2;, A, — a), \l/+(x, 2, A, a) = \l/_(x, z, A, — a), 
K+{dJ,\) = K4-a, -/3, A), 

while the t-evolution of K± is given by: 

K±{a, 0, A, t) = K±{a, 0, A, 0)e'^'^"^~'^^^'. (65) 

6 Dressing 

Dressing schemes corresponding to the three inverse problems presented in 
this paper can always be formulated (see, f.i., [12] for a general treatment). 
Here we restrict our attention to that corresponding to the RH inverse prob- 
lem of Section 5.1, but for the eigenfunctions 7i± defined in (j43|) . The proof, 
quite standard in the dressing philosophy, is left to the reader. 

Consider the following nonlinear RH problem on the real X-axis: 

7f+ = 7f_ + ^(7f_,A) =7l(ff_,A), AgM, (66) 

for the functions 7[+{x, z,t, X) and 7r_{x, z,t, X), analytic respectively in the 
upper and lower halves of the complex plane X, with asymptotics: 

- /- XN f X - Xz \ Q-i-{x,Z,t) „,,_2n IM -, /n^\ 

7r±(a;,^,A)= (^ ^_^^ J+^^A_ L + o{X^). |A| » 1, (67) 

where the vector R{n_, A) depends on {x, z,t) only through n_. Assume also 
the unique solvability of the RH problem / I6'6|) and of its linearized version 



z7+ = /7_ + p(7f_, A)z7_, A G M, 
(pk(C,A) = ^ 



dRdc,x) (68) 



Then tt-i- are solutions of Litt± = L2n± = 0, where Li, L2 are defined in (|i^| ) 
with: 

ui{x,z,t) = -Q±^, U2{x,z,t) = -Q±y. (69) 
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It follows that the potentials Ui^2 ore solutions of the nonlinear system of 
PDEs M'J\) . If, in addition, the spectral data satisfy the Hamiltonian con- 
straint 

Rk, + R2C, + {Ri,R2}^=0, (70) 

then the potentials Ui^2 satisfy the heavenly equation (QJj. At last, the reality 
constraints 

7l(7l(C',A),A) = C, VC, AeK (71) 

imply that 7?+ = ??_, and the reality of the potentials: Ui,U2 G M^. 
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